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ON THE EVALUATION AT {-i,i) OF THE TUTTE 
POLYNOMIAL OF A BINARY MATROID 

R.A.PENDAVINGH 



Abstract. Vertigan has shown that if M is a binary matroid, then |Ta/ (— i, i) \ , 
—i ^ the modulus of the Tutte polynomial of M as evaluated in (— i, i), can be ex- 

pressed in terms of the bicycle dimension of M. In this paper, we exactly 
determine TM{—i,i), and show how to evaluate this number in polynomial 
time. In particular, we describe how the argument of the complex number 
7a/ (—*!*) depends on a certain Z/4Z- valued quadratic form that is canoni- 
cally associated with M. 



o 
u 

i-H ■ 1- Introduction 

The Tutte polynomial of a matroid AI — {E,T) with rank function r is 

(1) rM(x,y):= ^(x~ir(^)-(^)(y-l)l^|-'-(^) 



FCE 



^ ■ Extending the work of Jaeger, Vertigan and Welsh [6], Vertigan investigates the 

^— ^ I complexity of r°(A^F,a;,y), the problem of evaluating TM{x,y) given x,y and a 

Qv ■ matroid M that is linearly represented over F [12]. Let j := e*^ denote a complex 

f^ , third root. 

CO . Theorem 1 (Vertigan, 1998). Let¥ be a finite field and {x,y) € A^ he a pair other 



than (0,0), (1,1), and such that [x — l){y -- 1) y^ 1. Then the problem t^ {A4r , x , y) 
is #P- complete, except when 

• |F| = 2, and (x, y) is one of (-1, -1), (0, -1), (-1, 0), (i, -i), (-i, i); 

• |F| = 3, and {x,y) is one of {j,f), {f,j); 

^! • \W\^4:, and (x, y) is (-1,-1) 

H ■ 

Cu ■ In the present paper we derive, for a binary matroid M, an explicit expression 

for Tm(— i,i) that can be evaluated in polynomial time. 

To the best of our knowledge, the complexity status of r°(A^F, 1, 1), i-e. counting 

the number of bases of an F-linear matroid, is open for |F| = 2,3,4. Otherwise, 

the above theorem is now complemented by explicit, polynomial-time computable 

expressions for the value of the Tutte polynomial on each of the special points: 

• for any matroid M, TM{x,y) is trivial to compute if {x,y) — (0,0) or if 
{x - l){y - I) = I. 

• for binary matroids M, |Tm(— i, i)\ was determined by Vertigan [12]. Below, 
we derive an explicit expression for TM{—i,i)- 

• for ternary matroids M, Jaeger [7] has determined \Tm{J,j'^)\, and Gioan 
and Las Vergnas 2 found Tm{J,P)- 
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2 R.A.PENDAVINGH 

• for quaternary matroids M, Tm{—^,—^) was found by Vertigan [12j . ex- 
tending a result for graphic/binary matroids by Rosenstiehl and Read [10] . 

The original motivation for this research was a computational problem that arose 
when writing a matroid package for Sage [11]. Testing whether two matroids are 
isomorphic can be made more efficient in practice by comparing matroid attributes 
that are invariants under isomorphism, avoiding more involved computation if the 
values of the invariants do not match. For general matroids, there are few such 
invariants that are polynomial-time computable. For binary, ternary and quater- 
nary matroids however, the values of the Tutte polynomial in the above-mentioned 
special points are clearly the kind of isomorphism invariant we can use for this 
purpose. With this objective in mind, we shall argue that computing TM{—i,i) as 
well as several related invariants takes 0{r{M)'^\E\) time all together. 

2. Preliminaries 

2.1. Matroids. A matroid is a pair M = {E,I) such that E is a. finite set and I 
is a set of subsets of E satisfying 

• 0eX 

• if / e I and J C /, then J e I 

• if /, J e I and |/| < | J|, then there exists an e e J \ / so that / U {e} e I 

A set i^ C i<^ is an independent set oi M ii F (£ X, and F is dependent otherwise. 
The rank oi a set F C E is r(F) :— max{|/| | / C F, / e I}. We generally follow 
Oxley [8] in our matroid notation, with the exception of the definition of linear 
matroid below. 

2.2. Linear matroids. If _E is a finite set, F is any field and F is a linear subspace 
of F^, then V determines a matroid M{V) := {E,I{V)), where 

(2) I{V) -.^ {F CE \ there is no u G F such that supp(w) C F} 

Here supp(w) := {e e E \ v,, ^ 0}. Thus r{M{V)) = \E\ - dim(\/). We will denote 

(3) V/e := {v\E-e \veV}, 

so that M{V/e) = M{V)/e. A matroid M = {E,I) is Imear over F if Af = M{V) 
for some linear subspace V C F-^; then V represents M. We write Aif for the set 
of matroids linear over the field F. Binary, ternary and quaternary matroids are 
matroids linear over GF{2),GF{3) and GF{'i), respectively. 

2.3. The bicycle dimension. Let y be a linear subspace of F-^, where |F| is one 
of 2, 3, 4. The bicycle dimension of V is 

(4) diV):^dim{VnV^). 

Here V^ := {w E ¥^ \< w,v >— for all v E V} as usual, where we take 
< w,v >:= '}2i'WiV, if F = GF{2),GF{3) and < w,v >:= J2i^i^'.^ where * : 
GF{A) — >■ GF{4:) is the unique nontrivial field epimorphism. 

Vertigan shows that if F be a linear subspace of F^ and |F| is one of 2, 3, 4, then 
d{V) depends only on M{V). 
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2.4. Quadratic forms. Let F be a finite-dimensional linear space over a field F. 
Let b : V xV ^>- ¥ he a, bilinear form. Then 5 : y — > F is a quadratic form associated 
with b if 

(5) q{Xx + fj.y) = X'^q{x) + fJ.'^q{y) + A/i6(x, y) 

for all A, /z G F and x,y eV. 

We say that a quadratic form is singular if the associated bilinear form is degen- 
erate, i.e. there exists a nonzero w & V such that b{v, w) ~ for all v ^ V . A 
basis ui , . . . , life of y is orthogonal with respect to b if 

(6) b{vi,Vj) = iii^ j 
The following is well-known. 

Lemma 1. If b : V x V ^ ¥ is a nondegenerate bilinear form and the field F has 
characteristic other than 2, then there exists a basis of V that is orthogonal with 
respect to b. 

Two quadratic forms q, q' on V are isomorphic if q{v) — q'{L{v)) for some linear 
bijection L : V -^ V. li vi, . . . ,Vk and wi, . . . ,Wk are 6-orthogonal bases of V and 
q is a quadratic form associated with 6, then Yli livi) is a quadratic residue if and 
only if Yli Q{wi) is a quadratic residue. Let xio) G {+I7 ^1} indicate whether the 
product is a quadratic residue. 

Theorem 2. If q,q' are nonsingular quadratic forms on V over a field of charac- 
teristic other than 2, then q is isomorphic to q' if and only if xil) ~ x{q')- 

The case when the characteristic is 2 is somewhat more involved. A basis 
vi, ■ ■ ■ , V2m of V is alternating with respect to b if 

,_s ,, \ f 1 if |«, ?'} = {A:, fc -I- 7n| for some fc e {1, . . . ,?7i} 

(^) ^(---.)-| otherwise 

Lemma 2. If b :V xV ^^¥ is a nondegenerate bilinear form and the field F has 
characteristic 2, then exactly one of the following holds: 

• V has a basis that is orthogonal with respect to b, or 

• V has a basis that is alternating with respect to b. 

Brown generalizes quadratic forms over GF(2) to Z/Alj-valued quadratic forms 
q : V 1-^ Z/4Z, satisfying 

(8) q{x + y)^q{x)+q{y) + a{b{x,y)) 

where b : V x V y~> GF{2) is a bihnear mapping, and a : GF{2) — >• Z/4Z is the 
unique nontrivial ring homomorphism (so a(l) = 2). Such a quadratic form q is 
nondegenerate resp. alternating if and only if the associated bilinear function b is. 
Brown (see [13| ) defines an invariant cr[q) such that 

(9) Y. *'^'^ = ^ 

x&V 

Wood [T^ has classified the Z/4Z-quadratic forms as follows: 

Theorem 3 (Wood, 1993). Ifq,q' are nonsingular "L/Alj-valued quadratic forms on 
V , then q is isomorphic to q' if and only if 

• (j{q) = cr(q'), and 

• q is alternating <^ q' is alternating. 



-dim(V) xig(g) 

e 4 
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3. A SPECIAL POINT OF THE TUTTE POLYNOMIAL 

3.1. Characterization of TM{—i, «)• Let V^ be a linear subspace of GF{2)^ , and 
let qv -.V ^ Z/4Z be defined by 

(10) qv{x) ■— |supp(x)| mod 4 

for alia; € V. Then qv{x + y) = qv (x) + qv (y) + a{b{x , y)) taking b{x,y) = Y.iXiyi, 
so that qv is a Z/4Z-valued quadratic form on V. 

If 1v{y) = for all 2/ e F n V^^, then for all x e F and all y G F n F-'- , we have 

(11) qv{x + y) ^ qv{x) + qv{y) + a{b{x, y)) = 17^(2;) 
and so then, we may define qv : V/{V n V^) -^ GF{2) by setting 

(12) qv{x + VnV^)^qv{x). 
Then q is nonsingular by construction. 

Theorem 4. Let V be a linear subspace of GF{2)^ and let M := M{V). Then 

iiV) 



(13) Tm{~i, i) = eT(-(9»+l^h3r(iS))y2" 
*/ 1v{x) = for all X (E V O V-^ , and Tnji—i, i) — otherwise. 

Proof. By an application of Greene's formula [4 (as in 12, p. 390]), we have 

(14) ^i«-(^) =z'-(^)(l~i)l^l-'-(^)TM(-i,*) 
xev 

Rewriting, we obtain 

(15) TMi-^, i) = V2"""+''<'''e^(l^l-3'-(s)) ^ ,«v(.) 

x£V 

If qy{x) = 2 for some xeVC^V^, then Y.xev *'"'^'''' = and hence TA./(-i, i) = 0. 
If on the other hand qv [x) = for all x e F n V^ , then by ([121 we have 

(16) \^ i<iv(^) ^ y" y^ j9v(i^+«)) ^ y^ y^jgX«+vny^) 

where V is any subspace of V so that F = (V^ n V^) F. The summation over 
V n F"*^ amounts to a factor 2''(^\ and using (jH]) on the non-degenerate form qv 
we obtain 

(17) Y. *''''^ = 2'^(^)y2'''"^^^e^-(^-) 

Substituting this expression in ([15]), and using that \E\ — r(Af ) = dim(y) = rf(V^) + 
dim(l/), we obtain ([HI). □ 

For comparison, we state the characterization of r(j, j^) due to Gioan and Las 
Vergnas j2 in similar terms. For a subspace V C GF{3)^ , let gy : T^ — ?► 6*^(3) be 
defined by 

(18) gy : a; H> |supp(x)| mod 3 

and let qv : V/{V n l^-^) ^ GF{3) be defined by qy(a; + Vr\V^)^ qv{x). 
Theorem 5. Let V C GF{3)^ be a linear subspace and M :— M{V). Then 

(19) TMij,f) = (-l)^^^^/l^l-'^(*^)(^y3)'^(^) 
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3.2. Complexity of TM{-i,i)- In what follows, let V C GF{2)^ be a linear 
subspace of dimension k. A q-basis is a basis vi, . . . ,Vk oi V such that 

• vi, . . . ,Wfc-(i(v) is an orthogonal or alternating basis of some subspace V 
such that V ^VQiVnV-^) 

• Vk-d{v)+i, ■ ■ ■ ,Vk is a basis of F n V-^ 

• qv(vk-d{v)+i) = 0, . . . , qv{vk-i) = 
Standard linear algebra techniques yield: 

Lemma 3. Given any basis ofV, computing a q-basis takes 0(dim(y)^|_E|) time. 

The following is straightforward from the definition of (/-basis and 0. 

Lemma 4. Let vi, . . . ,Vk be a q-basis of V. Then qv{v) = for all v £ V O V^ if 
and only if d{V) ~ or qv{vk) = 0. If so, then 

• if Vi, . . . , Vj._My\ is orthogonal, then 

(20) (T{qv) = W\qv{v^) = l}-^{i\qv{v^) = ?,} mods 

• if vi, . . . ,Vk-d(v) *s alternating, then 

(21) a{qv) = 4:#{ie {l,...,m}\qv{vi) = qv{vi+m)} mod 8 

, k-d(V) 

where m = — j — - 

Theorem 6. Let V C GF(2)^ be a linear subspace. Given any basis of V , the 
evaluation of Tj^fiy\{—i,i) iafces 0(dim(V^)^|i?|) time. 

Proof. To compute T'm(i/)(— i,i), it suffices to determine |_E|, rank(M(F)) — \E\ — 
dim(y), whether qv{v) = for all w e V^ n V^ , and if so, cr{qv)- Given a q-basis of 
V, this takes in 0{Aim{V)\E\) time. D 

As dim(y) = \E\ - r{M{V)) = r*{M{V)), this amounts to a complexity bound 
of 0{r*{M{V)f\E\) for evaluating TM(y)(-i, J) from a basis of V. We note that 
as in general Tm(x, y) — Tm* (y, x) = Tm* (y, 2^), we have 



(22) TM{v){-i,i) = TM(v^)i-h'i-) 

We may determine the latter in 0(dini(y-'-)^|i?|) = 0{r{M{V))^\E\) time from any 
basis of V^ . 

3.3. The canonical tripartition. Let V C GF{2)^ . We consider 

(23) F,:={eeE\ d{V/e) ^ d{V) + i}. 
Lemma 5. Let vi, . . . ,Vk be a q-basis of V . Then 

• F-i — [jLk-d{v)+i suppiv,) 

• Fi := supp(^-~-^^ '^i)\ P-i */ ^^6 q-basis is orthogonal, Fi — % otherwise 

• Fo = S\(F_iUFi) 

In particular, Fi = $ for any i ^ { — 1, 0, 1}. 

Proof. Let A be any k x E matrix over GF{2) such that V = rowspace(^). Then 
IZ-L = kernel(A), and d(y) = dim(y n V^) = k - rank(ylyl'^). If we write 

(24) A" := A[{1, ...,k},E-e] and a" := A[{1, ...,k},e] 

for e e E, then V/e = rowspace(^^) and d{V/e) = fc - rank(y4''(A'')^). Thus 

(25) d{V/e) - d{V) = rank(^A^) - rank(ylA^ + a''(a")^) 



6 R.A.PENDAVINGH 

Now consider the matrix A whose rows are the given g-basis of V. 

If e g supp(wi) for some i > k — d{V), so af ^ 0, then consider the matrix B 
that arises by adding the i-th row of A to the rows j G supp(a'^) — {i}. Then the 
rows of B again form a q-basis of V, b'^ is a unit vector, and hence d{V/e) — d{V) = 
v&nk{BB'^) - rank(SB^ + b^b'')'^) = -1. 

If of = for aU i > fc - d{V), then d{V/e) - d{V) = rank(^^^) - rank(ylA^ + 
a^a"")^) = rank(SB^) - rank(BB^ + b^b"")^) where B = A[{1, . . . ,k~ d{V)}, E]. 
So without loss of generahty, we may assume d{V) = 0. 

In case the rows of A are orthogonal, let B := v4[supp(a'^), £']. Then d{V/e) — 
d{V) = rank(AA^)-rank(AA'^ + a^(a'=)^) = rank(BB^)-rank(SB'^ + 6'=(fe^)'^) = 
rank(/) — rank(/ + J), where J is an all-one matrix with |supp(a'^)| rows/columns. 
Then d{V/e) - d{V) ^liiJ2^a■ = |supp(a^)| mod 2 == 1, and otherwise. 

In case the rows of A are alternating, it remains to show that d{V/e) — d{V) = 0, 
i.e. that AA'^ + a'^{a'^)'^ is nonsingular. If not, there is a nonzero vector x e GF{2)^ 
so that {AA^ + a^{a^Y')x = 0, or equivalently, so that AA^x = a^{a'^)'^x. As 
AA^x 7^ 0, we must have {a^Y^x = 1 and hence AA^x = a^ . As 



(26) AA' X 



/ 
/ 



{^kl2+\, ■ ■ ■ .'Xk-.X-y., ■ . ■ ,X^I2)^ 



this implies that ((f^x = 0, a contradiction. D 

The triple (F^\^Fq^F\) is known as the canonical tripartition of V. Originally 
this notion was developed for graphs by Rosenstiehl and Read [10], but the gen- 
eralization to binary spaces is straightforward. As is clear from the lemma, the 
determination of the canonical tripartition takes 0(dim(T^)|_E|) time given a q-basis 
oiV. 

3.4. A projection. Following Godsil and Royle[3], we call a subspace V C GF{2)^ 
pedestrian if the 'bicycle space' V n V-^ is trivial, i.e. if d{V) = 0. 

Pedestrian spaces are not rare. Fixing a subspace W C GF{2)^ of dimension k, 
there are exactly 2'^ distinct subspaces V C GF{2)-^~^'^ such that dim(V^) — k and 
V/e = W. Of these spaces V, exactly 2*= - 2^''^^^^ satify d{V) = d{W) - 1. From 
the remaining 2'=-'^(^), exactly half has d{V) = d{W), half has d{V) = d{W) + l. A 
straightforward analysis shows that ii \E\ tends to infinity, slightly less than 3/7-th 
of the subspaces V C GF{2)^ of dimension k will be pedestrian. 

If y C GF{2)^ is a pedestrian subspace then GF{2)^ = F ® F^, i.e any vector 
X £ GF{2)^ can be uniquely written as x — v + w, where v G V and w £ V^. 
Hence there is a unique linear map Try : GF{2)^ — > y so that Try (a;) G V and 
X — Try(a;) G V-^ for all x G GF{2)^. The matrix Qy so that 7ry(a::) — Qvx is 
textbook material in linear algebra; it is 

(27) Qv^A^{AA^)-^A 

where A is any matrix with independent rows such that rowspace(A) = V . That 
AA"'" is invertible follows from our assumption that V is pedestrian, that Qy G 
rowspace(^) = V is, clear, and we have x — Qyx G kernel(A) ~ V^ since 

(28) A(x - Qvx) = A{x - A^{AA^)-^Ax) ^ Ax - AA^{AA^y^Ax = 
It is not difficult to determine Qv from a g-basis of V . 
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Lemma 6. Let V C GF{2)^ be a pedestrian linear subspace, and let vi, . . . ,Vk be 
a q-basis of V . Then 

k fc/2 

(29) Qy = ^ v,vj orQv^Yl VtvI+k/2 + '"i+k/2vj 

i=l i=l 

if the basis is orthogonal or alternating, respectively. 
Proof. Let A be the matrix whose rows are vi, . . . ,Vk. Then 

/ 



(30) AA' = / or AA' = 



I 



depending on whether the g-basis is orthogonal or alternating. The result now 
follows from (1271). □ 



Let Gv be the support graph of the E x E matrix Qv, so V(Gv) — E and 

(31) E{Gv) = {ef I {Qv)ef ^ 0} 

Theorem 7. IfVC GF{2)^ and V' C GF{2)^ are pedestrian binary spaces, then 
M(V) is isomorphic to M{V') if and only if Gy is isomorphic to Gy . 

Binary matroid isomorphism may thus be reduced to graph isomorphism in the 
time it takes to write Gy, which is 0(dim(y)|ii^p). 

If we do not assume that V is pedestrian, then we may still write 

(32) (y + v^)/{v n v^) = v/{v n v-^) ® v-^/(y n v-^) 

Hence there is a unique linear map ny : {V + V^)/{V n V-^) -^ V/{V n V'^) such 
that 

(33) %(a;) G V/iVnV-^) anda;-^y(a;) e V-^/iVnV-^) 

for all X E {V + V-^)/{V D V^). Clearly, ttv coincides with Try if V is pedestrian. 

Lemma 7. Let V C GF{2)^ be a linear space and let V be any subspace ofV such 
that V = V +{Vn V^). Then nv{x + (V n V^)) = TTy{x) + [V n V^) for any 

xeV + V^. 

Proof Let x e V + V-^. li W C V-^ is such that V-^ ^ W ® (V n V^), then 
V + V^ — V <3)W <3) {V n V^) and we may write x = v + w + z, where v G V, 
weW, zeVnV^. Then 

(34) nvix + {V + V^)) = v + {Vr]V^) = 7r,>(x) + {V D V^) 

as required. D 

And so we have 

(35) nv{x+{VnV^)) = QyX+{VnV^) 

for any V such that V — V + {V O V-^). Such a V^ is determined as the span of the 
first k — d{V) vectors of a g-basis of V. 

If (-F-i, -Fo, -Fi) is the canonical tripartition of V , then from Lemma [5] we have 

(36) F_i = |J{supp(i;) I u e y n V^} 

Then (Q(/)e/ is independent of the choice V ii e, f E F := E\ F_i, since then 

(37) < ee + (y n V^),7rv{ef + {V D V^)) >= e^QyCf = (g^)e/ 
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Let Gv denote the support graph of Qy[F, F]. Then Gv — Gv if V is pedestrian. 
Theorem [7] thus extends to non-pedestrian spaces in a weaker form; 

Theorem 8. IfV Q GF{2)^ and V C GF{2)^' are binary spaces, then M{V) is 
isomorphic to M{V') only if Gv is isomorphic to Gy. 

4. Some remarks and conjectures 

4.1. Counting bases. As was mentioned in the introduction, the complexity of 
the foUowing problem appears to be open for any fixed field F: 

given: a basis of a subspace V C F^ 
find: the number of bases of M{V) 

Vertigan announced the resolution of this issue in J12j . but the paper he refers 
to remains unpublished. It seems natural to conjecture that this problem is #P- 
complete. 

In contrast, the matrix-tree Lemma states that if M is a regular matroid repre- 
sented by a totally unimodular matrix A with independent rows, then the number 
of bases of AI equals det{AA^). This result generalizes to sixth-root-of- unity ma- 
troids and even quaternionic-unimodular matroids (see [9] for the extension of the 
matrix-tree lemma to quaternionic matroids, as well as a noteworthy use of the 
matrix Qy in relation to counting bases in minors of M{V)). 

A bold conjecture one might phrase is that any minor-closed class of matroids 
of quadratic growth rate (see [1]) is such that the number of bases of a matroid 
in the class is computable in polynomial time from some succint description of 
the matroid. As leaving the nature of the succinct description open make this 
conjecture a bit too vague, we will be more specific. We conjecture: if P is a partial 
field so that the class of matroids representable over P has quadratic growth rate, 
the the problem 

given: a P-matrix A 

find: the number of bases of M{A) 

is polynomial-time computable. In particular, this conjecture would imply that 
counting the bases of a dyadic matroid takes polynomial time. 

4.2. Isomorphism testing for binary matroids. We consider the problem: 

given: bases for subspaces V C GF{2)^, V C GF{2)^' 
decide: if M{V) is isomorphic to M{V') 

This problem properly contains isomorphism testing for 3-connected graphs, and by 
a simple reduction the general graph isomorphism problem. The complexity of the 
latter problem remains open to this day. In practice one's best option for solving 
graph isomorphism is to use a canonical labelling algorithm as was described and 
implemented by McKay. 

We described several isomorphism invariants that one could compute from a 
basis of V in just 0(dim(F)|-Ep) time: 

• Tmcv) (-■*,«) 

• the cardinalities of the Fi in the canonical tripartition (i^_i,i^o,^i) 

• the number of edges of Gy, G[Fq\ and G[Fi\ 
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In a random selection of 10000 subspaces of GF{2)^'^ of dimension 10, all but 324 
of the pairs of were revealed to generate non-isomorphic matroids by a comparison 
of these invariants, or were identified as isomorphic by an application of Theorem 
[T] This means that in less than 1 in 100.000 cases, it was necessary to revert to 
other methods for testing isomorphism. As the time to compute the invariants is 
comparable in practice to the time it takes to compute the rank of a few sets X in 
M(V) from a basis of V^ this results in a sharp reduction in the average time taken 
by the isomorphism test. 

In a forthcoming Sage package for matroid computation, this technique has been 
implemented to speed up the isomorphism test for binary matroids, and similar 
methods have been implemented for ternary and quaternary matroids. 

Haggard, Pearce and Royle ^ describe a practical algorithm to compute the 
Tutte polynomial of a graph, which makes extensive use of graph isomorphism 
testing to reduce the overall computational effort. A possible application of our 
isomorphism test would be the extension of their method to a practical algorithm 
for computing the Tutte polynomial of a binary matroid. 
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